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We show how a quantum walk can be implemented for the first time in a quantum quincunx 
created via superconducting circuit quantum electrodynamics (QED), and how interpolation from 
quantum to random walk is implemented by controllable decoherence using a two resonator system. 
Direct control over the coin qubit is difficult to achieve in either cavity or circuit QED, but we show 
that a Hadamard coin flip can be effected via direct driving of the cavity, with the result that the 
walker jumps between circles in phase space but still exhibits quantum walk behavior over 15 steps. 

PACS numbers: 03.67.Ac, 42.50.Pq, 74.50.+r 



I. INTRODUCTION 

The quantum walk (QW) is important in physics as 
a generalization of the ubiquitous random walk (RW), 
which underpins diffusion and Brownian motion. QWs 
are important in quantum algorithm research [l[ because 
they exponentially speed up hitting time in glued tree 
graphs Q . Although realization of a QW by a quantum 
quincunx, analogous to the quincunx (or Galton Board) 
for realizing the RW, has been proposed in ion traps [3[ 
and cavity quantum electrodynamics (QED) [J, U, the 
former requires cooling of ions to the center-of-mass mo- 
tional ground state, and the latter hitting the atom with 
pulses that do not drive the cavity at all: these obstacles 
prevent quantum walks from being realized under fore- 
seeable experimental conditions (a classical optical simu- 
lation of a quantum quincunx has been performed [ffbut 
cannot be a proper QW without complementarity |7j). 

Here we devise a quantum quincunx that can realize 
a QW in the laboratory for the first time by: (i) em- 
ploying the Jaynes-Cummings model 0] to generalize the 
Hadamard transformation for coin flipping by directly 
driving the cavity rather than the atom; (ii) developing 
a theory of QWs on many circles in phase space (PS) 
rather than on a single circle as a consequence of gen- 
eralizing the Hadamard transformation; (iii) optimizing 
the protocol by having the duration of the generalized 
Hadamard transformation depend on the time-dependent 
mean photon number in the cavity (with detailed theory 
to appear elsewhere @); (iv) implementing in a super- 
conducting circuit QED system E| with a two-level 
Cooper Pair Box (CPB) serving as the quantum coin 
and a coplanar transmission line resonator with a single 
mode as the quantum walker; (v) introducing a double 
resonator scheme that can control decoherence while si- 
multaneously enabling strong coupling between the CPB 
and the microwave field and permitting fast readout; and 
(vi) using the Holevo standard deviation as a measure of 
phase spreading and showing that the rate of spreading 
can be tuned by controllable decoherence to observe the 



quadratic enhancement of phase spreading for the QW 
vs RW. 

In our scheme the QW is executed with indirect flip- 
ping of the coin via directly driving the cavity and allows 
controllable decoherence over circles in PS. Because the 
walker is directly driven rather than the coin, photon 
number is no longer conserved, and the walker jumps 
between circles in PS; however, a signature of quantum 
walking on circles in PS is evident in both the time- 
dependent phase distribution of the walker as well as via 
direct homodyne measurements to obtain the quadrature 
phase (QP) distribution for the walker. This signature is 
scaling of the standard deviation a (which measures the 
walker's spreading) that is linear in time oc t for the QW, 
and whose power decreases with increasing decoherence 
until attaining the classical RW scaling oc \ft for full de- 
coherence. This controllable decoherence is achieved by 
introducing a second low-Q resonator to obtain fast read- 
out [13 . 

II. BACKGROUND 

To understand the QW in PS, it is helpful to first un- 
derstand the RW in PS. The walker is a mode of the 
resonator, hence is equivalent to a harmonic oscillator, 
which can be described by its position x and momen- 
tum p. If the oscillator's energy E ~ (uj 2 x 2 + p 2 )/2 for 
unit mass and frequency w, then the oscillator follows 
a periodic circular trajectory of radius \f~E~ in PS with 
physical oscillatory motion 

x(t) = VEcosuot. (1) 

This oscillator can be modified to execute a RW on a cir- 
cle in PS by periodically applying an impulse that causes 
it to rotate either clockwise along the circle in PS by an- 
gle A9 or counter-clockwise by the same amount, with 
the choice of ±A6 strictly random. If A9 = 2ir/d, d e N, 
then the walker always remains on a (perhaps rotating) 
lattice on the circle with angular lattice spacing A6. 
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We refer to the coordinate in PS (x,p) as the walker's 
'location' in PS, and the coin flip randomness that de- 
termine clockwise vs counter-clockwise angular steps im- 
plies that the walker's location is indeterminate hence 
described by a distribution P(x,p). 

In an ideal QW on a circle, the walker+coin state is a 
density operator p € B(H) for H = H^^Hc (with walker 
space spanned by d — 1 discrete phase states [!, H| , 

H w = span{|# m — 2mir/d)} (2) 

and coin space H c = span{|0), |1)}), where B(H) is the 
Banach space of bounded operators on H. The walker's 
phase distribution on the circle is 

P w (0) = (6\p w \6)/2ir, Pw = Tr c p, (3) 

with d equally spaced values of 6 = 9 m . However, here 
the walker is following a circular trajectory in PS so the 
QW's Hilbert space is H w = span{|?i);n G N}, with \n) 
a Fock state (h eigenstate). 

Alternatively the generalized position representa- 
tion {1^)^} can be used with \x) v an eigenstate of 
x cos <p + psunp, for x, p the canonical operators satis- 
fying [x,p] = i (h = 1). The QP distribution is 

P v (x) = v (x\p w \x) v (4) 

with ip the phase of a local oscillator. A convenient corre- 
spondence between the classical and quantum PS trajec- 
tory of the walker is provided by the Wigncr quasiprob- 
ability distribution 

W(x,p) = J°° ^ m {x - y/2\p w \x + y/2) (5) 

whose marginal distributions are P v (x). Our scheme for 
realizing the first experimental QW builds on the cavity 
QED quantum quincunx pj], which alternately applies a 
coin flip Hadamard gate 

ff = |+)<0| + |-)(l|, (6) 

for 

|±) = (|0) ± |1))/V2, (7) 

followed by a rotation of the walker's state in PS by ±A8 
with the sign depending on the coin state (with n = d^d 
for a = [lux + ip] / ^/2lj) : 

F = exp (ina z A8) . (8) 

The initial state of the walker is a coherent state 
a£l and 

(n\a = V^o) = ye-^of^/nl. (9) 

Fig. [1] depicts PS, including how d is chosen given an 
initial walker state |a) 

n + \ffi < d < 2Try/n. (10) 
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FIG. 1: PS diagram depicting circles of fixed radius y/n] 
for circle j. The circle of radius |a| depicts 1/e contours of 
the Wigner function for coherent states separated by angu- 
lar spacing 2n/d to make them distinct. A generic Wigner 
function 1/e contour (boundary of shaded areas) is depicted 
for circle j with mean (n) — fij. Solid boundaries represent 
contours around positive peaks, while dashed boundaries are 
for contours around negative peaks. 



With this choice, the walker's angular step size is large 
enough to ensure enough distinguishability to yield a cir- 
cular QW signature. 

The simplicity of this model is rooted in the commu- 
tativity of (FH) N with n, hence a constant of motion; 
thus the walker's distance from the PS origin \ffi is fixed. 
Unfortunately alternating between F and H is not prac- 
tical in circuit QED. For the QW to be realized, and also 
to have controlled dccohcrcncc, a time-dependent driv- 
ing field for the resonator is needed. As we see below, 
the sacrifice is that [FH, h] ^ 0, but for realistic circuit 
QED conditions, the QW on circles is (surprisingly) evi- 
dent provided that the step-by-step H pulse duration is 
judiciously chosen. 

The goal is to observe the spreading of the QW's phase 
distribution, and the signature of the QW is that this 
spread is linear in time: specifically the standard devi- 
ation for the phase distribution satisfies a oc t for the 
QW, whereas a oc Vt for the RW, and the QW can 
be tuned continuously from the RW by controlling de- 
coherence. As phase is periodic, the usual root-mcan- 
square approach to a is problematic; instead we employ 
the Holcvo standard deviation [l3| 

a H = v l<c^)|- 2 -l (11) 

for 

(e^) = / d(j>P{^ (12) 
Jo 
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with respect to any phase distribution P(<fi) (0). The 
Holevo standard deviation is equivalent to the root-mean- 
square definition for small spreading on the circle, and gp 
naturally quantifies dispersion over all <fi <E [0, 2tt) 14|. 
Thus, for sufficiently short times (less than the time that 
phase distribution spreads over a significant fraction of 
the circle in PS), if the relation between phase spreading 
on the circle with time is a power law, then 



rncr H = ?lni + £ (13) 
with ? = 1 for the QW and q = 1/2 for the RW. 

III. SUPERCONDUCTING CIRCUIT 
QUANTUM ELECTRODYNAMICS 

The circuit QED Hamiltonian is 

H = H JC + H d , H d = e(t) (aV^' + ae^) (14) 

for H d the time-dependent driving field Hamiltonian and 

Hjc = w," + ^a<7 z /2 + g{a) a~ + aa+) (15) 

the Jaynes-Cummings (JC) Hamiltonian with uo a and w r 
the coin and resonator frequencies, respectively, and g the 
coupling strength. It is sufficient to let e(t) be a square 
wave so e is a constant (e = when the field is off). In 
the dispersive regime, 



|A| = |u>„ - w r | > g, 



(16) 



and in a frame rotating at uj& for the qubit and the res- 
onator, H can be replaced by the effective Hamiltonian 

HcS =xna z / - 8 Aa a z /2 

- 5 dl h + Q, R a x /2 + e(a t + a) (17) 



with 



fl R = 2ge/S dr 
the Rabi frequency, and 

X = ,9 2 /A 



(18) 



(19) 



(20) 



the cavity pull of the resonator. 

The first term in the above expression effects the coin- 
induced walker phase shift. The atom transition is an 
ac-Stark shifted by g 2 hA. To implement 

H = exp [itH^R^/2] (21) 

on the coin, we choose 

uj d = 2ng 2 /A - 2ge/A + uj & (22) 



with pulse duration in = k/2£Ir and Or, a function of 
average photon number 



The free evolution 



n(t) = Tr(np w ). 



exp (-iH cS t) 



(23) 



(24) 



continues even when the driving field is off (e — 0). For 
time r between i7-pulses, the walker steps through an 
angle 



A6 « ±g\r + t H )/A. 



(25) 



Whereas the ideal QW conserves photon number, 
Eq. (jTTJ) violates this, which we interpret as the walker 
wandering between circles in PS. Circles have radii 
\J n(t), and Ud and £r adjusted due to n(t) to ensure 
that the angular step size Ad is constant regardless of 
how far the walker is from the PS origin. The mean pho- 
ton number n(t) can be calculated as follows [9(. We 
begin by solving the Schrodingcr equation 



dt 



= -iH c s \<f) 



from time t = to t = N(t^' + r), for N the number of 
steps 

4 0) = tt [A + 2(| a | 2 + 1) X - 2.96/ A] /Age (27) 



(26) 



and beginning with the initial state 

|^ o ) = (|0)+i|l))|a)/V2. 



(28) 



Figure 2 plots n{t) for a = 3, d = 21 and realistic 
system parameters 

(w a , u T , 5 , e)/27r = (7000, 5000, 100, 1000)MHz. (29) 

It is evident in the figure that the mean number oscillates 
and then settles down during the free evolution so the 
walker is concentrated on a circle of squared radius 



-i «(4 0) +-) 



0-i)(*^+t) 



n(t)dt (30) 



at step j. The corresponding Hadamard pulse dura- 

(1) 

tion fjj for each step j is 



(.;) 



n[A + 2(n j + l) X -2ge/A]/4ge (31) 



for j e N, n Q = 9 and t$ = 0.01567 /is. 

Alternatively, to compensate for photon number fluc- 
tuations, it is possible to vary the frequency of the 
Hadamard pulse rather than its duration. For large step 
number N, the photon number distribution is found nu- 
merically to approximate 



P(n) 



{n\ p w \n) 



- n n n /n\, 



(32) 



4 



and the width of the photon number distribution P(n) is 
closely approximated by 

5n = y/(h 2 } - (n) 2 w yffi (33) 

in the numerical simulations. Thus, number spreading 
is negligible, and the walker can be regarded as indeed 
being concentrated in the locality of one circle in PS [§]. 
As the walker's initial state is a coherent state with a 
Poissonian number distribution, the fact that the width 
remains Poissonian, and the amplitude is relatively con- 
stant, indicates that the initially well localized walker 
continues to be localized with respect to amplitude in 
the phase space over time. 

IV. LOCALIZATION OF WALKER IN PHASE 
SPACE 

We observe from numerical simulations that the 
walker's location in phase space is effectively localized to 
a circle in phase space with radial width given by (f3"3")) . 
Here we explain why this confinement to the vicinity of a 
circle in phase space with radius \ffi is reasonable. This 
discussion is based on a theoretical analysis of quantum 
walks on circles in phase space . 

The evolution of the joint walker+coin system is gov- 
erned by the effective Hamiltonian (jT7|) . There are five 
terms on the right-hand side of this Hamiltonian. Let us 
understand each term and its effect on the dynamics to 
see why the walker is localized to a circle in phase space. 



f . The first term involves na z , which is responsible for 
entangling the evolution of the coin and the walker, 
effectively to make the walker evolve clockwise or 
counterclockwise at constant amplitude with the 
orientation entangled with the state of the coin. 

2. The second and third terms involve the operators 
o z and n, respectively, which correspond to ener- 
gies, hence frequencies, for the coin and walker. 

3. The fourth term, involving <j x , is responsible for 
the Hadamard coin flip and is proportional to the 
Rabi frequency, which is itself proportional to the 
pulsed driving field e (jT9j) . 

4. The fifth term is a displacement involving the oper- 
ator a + d^, which pushes the walker off the circle, 
and is proportional to e. 

In making the quantum walk work, the goal is to make 
the Rabi frequency S1r large but keep e small in order 
to flip the coin but minimally shift the walker from the 
circle. 

We have been able to simultaneously achieve the two 
conditions of large Rabi frequency and small displace- 
ment. By meeting these two conditions, the evolution is 
closely approximated by the unitary evolution 

H g) D(i\/V2) (34) 
for the quantum walk on circles in phase space @ , with 
A = \[2(±Yi the size of the displacements from the circle. 
The resultant photon number spread after N steps is @ 



+ cos AO) 



2 + cos A6 - cos NA6 



cos(l-iV)Ae 
sin 2 A6/2 



+ cos AO 



(35) 



r 



For large no and small A, 5n ~ s/fi,. Hence, in the 
asymptotic large mean photon number n limit, the re- 
duced walker state has support almost entirely from co- 
herent states with amplitude \ffi. This means that the 
joint state of the walker+coin can be regarded approxi- 
mately as an entanglement of a walker in superpositions 
of coherent states with the coin state. This approxima- 
tion guarantees that the walker can be regarded as being 
localized to one circle in phase space with a Poissonian 
spread in photon number that does not increase signifi- 
cantly over time provided that the phase steps are small 
and n is large. 



V. OPEN SYSTEM AND MEASUREMENT 

Coupling to additional uncontrollable degrees of free- 
dom leads to energy relaxation and dephasing in the sys- 



tem. In the Born-Markov approximation, these effects 
can be characterized by a resonator photon leakage rate k 
(determined at fabrication time by the resonator input 
and output coupling capacitances), an energy relaxation 
rate 71, and a pure dephasing rate 74, for the qubit. The 
open system thus evolves according to 



P = -1 
with 



H e g,P 



KV[a]p + 7iX>[*_]/> + { lv /2)V[a z ]p, 

(36) 



v[L] P = (2L P D - Dip - pDi)/2. 



(37) 



Relaxation and dephas ing of a charge qubit in circuit 
QED were reported in [ll| as T\ = 7.3 ^s and T 2 = 500 
ns. These translate to 



7i/2tt = 0.02MHz 



(38) 
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FIG. 2: Average photon number n vs evolution time t (/is) 
with period + r). 

and 

l4> /2n = ( 72 - 7i/2)/2tt = 0.31MHz. (39) 

The master equation (|36j) is used to compute p(t) from 
which the reduced state of the walker p w is obtained and 
thence the phase distribution P w (9) ([3]). Empirically the 
phase distribution can be obtained by performing full 
optical homodyne tomography on the transmission line 
resonator to obtain W(x,p) [17| from which P w (8) can 
be computed and o"a(t) thereby determined. The scal- 
ing of cth with t in Eq. (|13| is a convenient empirical 
signature of the RW vs the QW. 

As full tomography is expensive, it would be valuable 
to observe a QW directly from optical homodyne detec- 
tion with a single choice of local oscillator phase p instead 
of having to scan over many ip for full tomography. Our 
simulations show both an(t) for P w (6) and ctqp for the 
QP distribution at tp = 0, with the choice of ip = cor- 
responding to having a local oscillator that is in phase 
with the walker at t = 0. For our choices of realistic 
experimental parameters, the walker's phase distribution 
spreads from — tt to it on a scale of 15 steps so simulations 
are limited to fewer than 15 steps before the spreading 
effectively saturates and our comparison of the QW vs 
the RW breaks down. 

Simulated evolution of <th(£) and (JQp(t) are presented 
in Fig. G3a,b), respectively. Corresponding linear re- 
gression data presented in Table U clearly reveals slopes 
compatible with the characteristic quadratic decrease in 
phase spreading for increasing decoherence of the QW 
until the transition to the RW 0]. 

These results show the significance of k in decoherence 
from the QW to the RW. Moreover k is much more im- 
portant than 71 and 7^ with respect to the scaling of 
ch,qp with t. The pure dephasing rate 7^ mainly leads 
to smearing of the phase distribution and the phase dis- 
tribution loses its symmetry. Furthermore, the effect of 
the energy relaxation rate 71 is small here compared with 
k because of our realistic choice of parameters. 

Unfortunately k must be low to obtain a QW yet high 
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ln(t)[ln(jis)] 

FIG. 3: (Color online) Ln-ln plot of walker spread a vs t 
(us) for (a) phase distribution on the circle in PS and (b) 
QP distribution, for different k and fixed (0,71, 74,) /2n = 
(100,0.02,0.31) MHz. 
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ln<r H 


Alncr H 


r 





0.924 


0.009 


0.442 


0.004 


0.990 


0.05 


0.879 


0.013 


0.362 


0.006 


0.991 


0.1 


0.822 


0.014 


0.279 


0.008 


0.992 


0.3 


0.615 


0.022 


-0.025 


0.012 


0.993 


0.5 


0.447 


0.030 


-0.309 


0.017 


0.990 



TABLE I: The linear regression data In cth = (s ± As) hit + 
(In CTjj ± AhiCT H ) of the Holevo standard deviation of phase 
distribution in ln-ln scale for lossy cavities. 



to allow fast readout. This dilemma is resolved by in- 
stead using two modes [lH: one resonator (labeled a) of 
high-Q and acting as the walker and a second resonator 
(labeled b) of low-Q and used for fast readout [12| • The 
microwave radiation from resonator a is coupled into res- 
onator b, and measurements ensue on resonator b. Mea- 
surements must be quick on the time scale of the walker's 
steps, so it cannot be longer than the time scale between 
pairs of Hadamard pulses. In the two-resonator system, 
«b bas therefore a lower bound of 0(g 2 /A). 

Due to coupling to resonator b, the transition fre- 
quency of the CPB in the resonator a and the pure de- 
phasing rate are changing with the cavity field in res- 
onator b, so the master equation for resonator a and the 
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As 


In ctqp 


Alno^p 


r 





0.937 


0.006 


0.093 


0.004 


0.988 


0.05 


0.892 


0.009 


0.016 


0.006 


0.989 


0.1 


0.832 


0.013 


-0.068 


0.008 


0.990 


0.3 


0.634 


0.020 


-0.372 


0.012 


0.993 


0.5 


0.453 


0.034 


-0.677 


0.019 


0.990 



TABLE II: The linear regression data In ctqp = (s± As) lnt + 
(lno-Qp ± AlnoQp) of the standard deviation of QP distribu- 
tion in ln-ln scale for lossy cavities. 

CPB is modified to [T2] 

9 = -i[H a , p] + n & V[a]p + jiV[a.]p + l*±£^V[a z ]p, 

(40) 

with 

H s = w lh n+ ^jp-a z + [gtfa- + e^aV^* + he] (41) 
where 



In circuit QED, the amplifier thermal noise is signifi- 
cant, with more thermal photons present than the mean 
resonator photon number [ill ]. However the signal QP 
distribution and full Wigner function can be obtained 
from the resultant homodyne detection statistics by con- 
volving the readout with a thermal function (which is a 
Gaussian mixture of Gaussian states centered at the PS 
origin). The spread of the convolution is determined by 
the mean thermal photon number, which is typically 20. 

The result is expected to be noisy quadrature phase 
readout, but repetition will yield, on average, the de- 
sired linearity of lntr vs Int. We can use a filter al- 
gorithm [l^| which takes the inversion formula for the 
measured Radon transform of the Wigner function with 
thermal noise and reconstruct the Wigner function of the 
corresponding noiseless signal. Hence we can achieve the 
noiseless phase distribution and QP distribution with the 
Wigner function of the noiseless signal. 



VI. CONCLUSIONS 



n'(t) = u> & + 2 Xb (K(i)| 2 + 1/2), (42) 

and 

r m = 8 X l\a h (t)\ 2 /K h . (43) 

In these expressions, Xb = .9b/( w a — w r t>) i s the cavity 
pull of the measurement resonator b and a,h{t) is the 
classical part of the measurement cavity field. Under 
these conditions, the master equation for the CPB and 
resonator b (|40|) is identical to Eq. (|36|) except for a pa- 
rameter change. It is interesting to note that, in this two- 
resonator case, qubit dephasing can be tuned by changing 
the number of photons injected in the read-out resonator 
(corresponding to measurement-induced dephasing [Tel]). 
As a result, it should be possible to observe the cross- 
over between the QW and the RW as a function of this 
tunable dephasing. 

In summary we have introduced the following proto- 
col to implement the QW on circles in PS using circuit 
QED. (i) Solve the Schrodinger equation from t = to 
t = N(t { ^ + t), to obtain the mean photon number n{t) 

from which the sequence {nj} and {tjj } are obtained, 
(ii) Prepare a high-Q resonator in its vacuum state (by 
simply cooling), (iii) When the vacuum state is prepared, 
inject a microwave field with a Gaussian pulse shape and 
temporal width Tq in order to prepare the high-Q res- 
onator in a coherent state |a). (iv) After t = 2Tq, im- 
plement the Hadamard pulse on the CPB by injecting a 
square pulse of frequency lu^ into resonator b over time 
scale t\p for step j. (v) Terminate the external driving 
of the resonator to allow free evolution over time scale 
r. (vi) Repeat steps (iv) and (v) N times, (vii) Perform 
full tomography on resonator b by performing homodyne 
measurement over many values of cj> and use standard in- 
version technique on the data. 



We have shown that the QW on circles in PS can be 
implemented via circuit QED with and without open sys- 
tems effects using realistic parameters and find that the 
signature of the QW is evident under those conditions. 
Moreover full tomography may not be required because 
direct homodyne measurements over few quadratures re- 
veal an unambiguous QW signature, and the RW can 
controllably emerge by tuning decoherence. Our scheme 
shows how a QW with just one walker can be imple- 
mented in a realistic system for the first time, and con- 
trollable decoherence can be performed thereby allowing 
continuous tunability between the quantum and classical 
regimes. 

Our scheme allows only a finite number of steps before 
the quadratic enhancement in phase spreading breaks 
down. As the phase step A9 must be strictly greater 
than l/5n, the number of steps N has an upper bound 
because of the desire to avoid wrap-around effects (the 
walker going around the circle), then N < 2ir/ AO. There- 
fore, N < 2n8n ?s 2tt^/E provides an upper bound on N. 

Hence the number of steps has an upper bound that 
is determined by n. If the mean number of photons in 
the resonator is increased, so is the allowed number of 
steps. Physically, however, the mean number of photons 
cannot be too large because the dispersive approximation 
that we exploits fails for large cavity photon number: this 
breakdown occurs for critical photon number [IcL I20I ] 



Eq. (|4"4"|) yields an upper bound of photon number in 
the resonator and limits the number of steps the walker 
can take and exhibit a quadratic enhancement of spread- 
ing. In our example, A = 2000 MHz and g/2ir = 100 
MHz so n cr it = 100. Here we have treated the case of just 
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nine photons, and seen strong evidence of a quadratic en- 
hancement of phase diffusion, so the quantum quincunx 
effectively works well below this critical photon number 
where the dispersive approximation breaks down. 

In our scheme, the pulse duration in is adjusted each 
time according to the predicted mean photon number, 
but precise control of the Hadamard pulse may be diffi- 
cult to achieve in experiments. 

In Fig. 4 we have shown a simulation of the cases with 
and without adjusting the duration of the Hadamard 
pulse sequence. Not adjusting the Hadamard pulse du- 
rations still yields a quadratic enhancement of the phase 
distribution due to the quantum walk, but it breaks down 
earlier. In our simulation, the breakdown occurs af- 
ter 10 steps rather than 15 for the adaptive pulse du- 
ration. For the first 10 steps, the numerically simu- 
lated standard deviation for QP distribution and the 
Holevo standard deviation in ln-ln scale are respectively 
shown to be approximately linear in hiN: Iiktqp = 
(0.939±0.007) In A+ (-2.090±0.005), the r coefficient is 
0.99, and lncr H = (0.890±0.006) In A+ (-1.563±0.003), 
the r coefficient is 0.96. 



FIG. 4: (Color online) Ln-ln plot of walker's phase spread a 
vs N (the number of the steps) to N = 15 with fixed pulse du- 
ration 0.0157/xs for (a) phase distribution on the circle in PS 
and (b) QP distribution for the classical random walk (solid 
line) and the QW (red dots) with a — 3, d — 21 and typical 
system parameters («», u T , g, e)/2n = (7000,5000,100,1000) 
MHz. In comparison, the blue dotted lines are for the QW 
with adaptive pulse durations, which breaks down after 15 
steps. 
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